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ABSTRACT 

We analyze the molecular H 2 emission and the stellar kinematics in a data cube of the nuclear 
region of M104, the Sombrero galaxy, obtained with NIFS on the Gemini-north telescope. After a 
careful subtraction of the stellar continuum, the only emission line we detected in the data cube was 
H 2 A21218. An analysis of this emission revealed the existence of a rotating molecular torus/disk, 
aproximately co-planar with a dusty structure detected by us in a previous work. We interpret these 
two structures as being associated with the same obscuring torus/disk. The kinematic maps provided 
by the Penalized Pixel Fitting method revealed that the stellar kinematics in the nuclear region of 
M104 appears to be the result of the superposition of a “cold” rotating disk and a “hot” bulge. Using 
a model of a thin eccentric disk, we reproduced the main properties of the maps of the stellar radial 
velocity and of the stellar velocity dispersion, specially within a distance of 0"2 from the kinematic 
axis (in regions at larger distances, the limitations of a model of a thin rotating disk become more 
visible). The general behavior of the /13 map, which is significantly noisier than the other maps, was 
also reproduced by our model (although the discrepancies, in this case, are considerably higher). With 
our model, we obtained a mass of (9.0 ± 2.0) x 10 8 M e for the supermassive black hole of M104, which 
is compatible, at ler or 2<r levels, with the values obtained by previous studies. 

Subject headings: galaxies: active — galaxies: nuclei galaxies: individual(M104) — galaxies: kine¬ 
matics and dynamics — Techniques: spectroscopic 


1. INTRODUCTION 

Due to its proximity and brightness, M104 (NGC 
4594), the Sombrero galaxy, has been observed in many 
spectral bands. It is an SA(s)a galaxy with a very in¬ 
clined (nearly edge-on) disk, at a distance of 9.2 Mpc. 
Its optical spectrum shows a number of low ionization 
emission lines and, based on this property, the nucleus 
of this object has been classified as a Low Ionization Nu¬ 
clear Emission-line Region (LINER - Heckman 1980). 

Radio observations have shown that M104 liasa com¬ 
pact nuclear source (< lpc). [Hummel et all (11984H . us¬ 
ing observations of the continuum at 2.6 and 20 cm ob¬ 
tained with the Very Large Array ( VLA), concl uded t hat 
t he ce ntral source has a flat spectrum and [Baiaia et all 
(tl988ll . using data at 1.49 GHz obtained with the VLA, 
detected a variable continuum in that source. 

An X-ray co m pact nuclear source was also observed. 
iPellegrini et al.l (120021 1 analyzed Beppo SAX and Chan¬ 
dra data and verified that the hard X-ray spectrum of 
the central source is well described by a power law. The 
authors obtained a luminosity of L x = 1.2 —2.3 x 10 4Q erg 
s~\ in the spectral range 2-10 keV, for that source. All 
these radio and X-ray properties are compatible with the 
e xistence of an ac tive galactic nucleus (AGN) in M104. 

iNicholson et all (11998H . using data obtained with the 
Hubble Space Telescope (HST ) and the ASCA, concluded 
that t he model of photoionizaton by a low luminosity 
AGN (jFerland fc NetzeJ 119831 : lHalnern fc Steineil Il983fl 
is the most adequate to explain the emission-line spec¬ 
trum of this galaxy. In iMenezes et al.1 (12013H , we ana¬ 
lyzed a data cube of the nuclear region of M104, obtained 
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with the Integral Field Unity (IFU) of the Gemini Multi¬ 
object Spectrograph of the Gemini-south telescope, and 
reported the discovery of collimation and scattering of 
the AGN emission. We also detected a broad component 
of the Ha emission line, which is another evidence for 
the existence of an AGN in M l04. Such broad compo - 
nent was previously reported bv lKormendv et all ( 1996 ). 
using HST data. 

iBendo et al.l (12006H , using infrared data obtained with 
the Spitzer and with the James Clerk Maxwell Tele¬ 
scope of the Submillimetre Common-User Bolometer Ar¬ 
ray, verified that the spectral energy distribution of the 
nucleus of this galaxy reveals the existence of hot dust. 
Besides that, the authors also concluded that the low flux 
at 160 /irn seems to indicate that the LINER emission is 
a consequence of the lack of cold gas to feed a possible 
more intense Seyfert activity. 

The mass of the supermassive black hole (SBH) in the 
nucleus of Ml0 4 was o ne of the first SBH masses to be 
determined. iKormendvl (1198811 analyzed data obtained 
with the Canada-France-Hawaii Telescope (CFHT ) and 
found evidence for a SBH with a mass of 10 9 M^. Thi s 
same value was also obtained bv lEmsellem et aLl(jl994ll . 
using axisymmetric Jeans models (a generalization of 
the Jeans formalism is described by Cappe ll ari 2008) 
and a dditional CFHT data, and by iKormendv et al.l 
(fl996ll. usi ng HST and higher resolution CFHT data. 
iMagorrian et all (11998H . using also axisymmetric Jeans 
models, HST photometry, and kinematics obtained 
from ground based observations, determined a mass of 
(6.4718'ia) x 10 8 M e for the SBH. Finally, Uardel et all 
((Mill . using high-resolution HST spectra, long-slit spec¬ 
tra obtained with GNIRS of the Gemini telescope, and 
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also integral-field kinematics obtained with SAURON, 
appli ed the Schwarz schild method of superposition of or¬ 
bits (|Schwarzsc hild 1979) and obtained a SBH mass of 
(6.6 ±0.4) x 1O 8 M 0 . 

In this paper, we present an analysis of a data cube of 
the nuclear region of M104, in the K band, obtained with 
the Near-infrared Integral Field Spectrograph (NIFS) on 
the Gemini-north telescope. The use of the adaptive op¬ 
tics (AO) during the observation resulted in a high spa¬ 
tial resolution for the data cube, which also has a com¬ 
bination of high signal-to-noise and spectral resolution. 
In order to improve the quality of the results, we treated 
and analyzed the data cube with a variety of techniques. 
Our goals are: to analyze the emission-line spectrum of 
the nuclear region of All 04 in the K band and to ana¬ 
lyze the stellar kinematics around the nucleus, with the 
additional purpose of obtaining an estimate for the mass 
of the nuclear SBH. 

This paper is organized in the following way: in Sec¬ 
tion 2, we describe the observations, the data reduction, 
and the data treatment. In Section 3, we analyze the 
emission-line spectrum of the data cube. In Section 4, 
we analyze the stellar kinematics and apply a dynamical 
modeling, in order to reproduce the observed kinematic 
maps and to obtain an estimate for the mass of the cen¬ 
tral SBH. In Section 5, we discuss our results and make 
comparisons with previous studies. Finally, we present 
our conclusions in Section 6. 

2. OBSERVATIONS, REDUCTION AND DATA TREATMENT 

The observations of the nuclear region of A1104 were 
made on 2011 April 26, (Programme: GN-2011A-Q-83) 
with NIFS on the Gemini-north telescope. This instru¬ 
ment has a field of view (FOV) of 3" x 3" and works 
together with the AO module ALTAIR, which allows a 
maximum spatial resolution of 0"1. The final products 
from observations made with this instrument are data 
cubes, containing two spatial dimensions and one spec¬ 
tral dimension. Five 640 second exposures of the nuclear 
region of M104 were made, in the K band, with a central 
wavelength of 2.2 /im. The final spectra had a coverage 
of 2.01 - 2.42 /jm and a resolution of R ~ 5155. 

The following calibration images were obtained during 
the observations: flat-field, dark flat-field (used for dark 
current subtraction), Ronchi-flat (used to compute spa¬ 
tial distortions and to perform a spatial rectification), arc 
lamp and sky field. The A0V standard star HIP 56901 
was also observed to be used for the telluric absorption 
removal and for the flux calibration. The data reduc¬ 
tion was made with the Gemini IRAF package and in¬ 
cluded: determination of the trim, sky subtraction, bad 
pixel correction, flat-field correction, spatial rectification, 
wavelength calibration, telluric absorption removal, flux 
calibration (considering that HIP 56901 has a blackbody 
spectrum), and data cube construction. Five data cubes 
of M104, with spatial pixels (spaxels) of 0"05 x 0"05, were 
obtained. 

We applied a correction of the differential atmospheric 
refraction (DAR) to all the obtained data cub es, using 
the p ractical approach described in iMenezes et al.l 
pi . Then, in order to combine the five corrected 
data cubes into one, we divided them in two groups. 
Numbering the data cubes as 1, 2, 3, 4, and 5 (following 
the order of the observations), we applied the following 


division: 

Group 1: 1, 2, and 5 
Group 2: 1, 3, and 4. 

A median of each group was calculated and, after 
that, an average of the two obtained data cubes was 
calculated, resulting in the combined data cube. 

In order to improve the visualization of the contours 
of the structures, we performed a spatial re-sampling to 
all the images of the combined data cube (conserving the 
surface flux), to obtain spaxels of 0"021. The improve¬ 
ment of the appearance of the images in the data cube 
is not only reason to apply the spatial re-sampling. We 
verified that the final spatial resolution is higher when 
the Richardson-Lucy deconvolution (see more details be¬ 
low) is applied to a spatially re-sampled data cube. After 
that, we applied a Butterwortli spatial filtering to the re¬ 
sampled images of the data cube, in order to rem ove high 
spatia l frequency noise. As explained in Menezes ct al. 
(l20H . we used a Butterworth filter corresponding to the 
product between an elliptic filter and a rectangular one. 
The order of each one of these filters was n = 2 and their 
cutoff frequency was 0.19 Ny (Nyquist frequency). 

The data cube of M104 revealed the existence of an 
“instrumental fingerprint”, in the form of large verti¬ 
cal stripes across the images. This feature also has a 
specific low-frequency spectral signature. The cause of 
the instrumental fingerprin t on NIFS data cubes is still 
unknown. As explained in IMenezes et al. (12014! ), NIFS 
reads out through four amplifiers in 512 pixel blocks in 
the spectral direction. Therefore, a slight bias level dif¬ 
ference between these amplifiers could be an explanation 
for this fingerprint. We used the Principa l Co mponent 
Analysis Tomography technique dSteiner et al.l 12009 1 to 
remove the fingerprint fro m the data cube of M104, fol- 
lowing the prescription of lMenezes et al.1 (1201412015D . 

The last step of our data treat ment procedure was the 
Richardson-Lucy deconvolution (Richardson 1972; iLucvl 
1974), which was applied to all the images of the data 
cube, in order to improve the spatial resolution of the 
observation. The point-spread function (PSF) used for 
this process was an image of the data cube of the star 
HIP 56901. The use of an image of the standard star 
data cube was necessary, because we were not able to 
determine any other reliable estimate for the PSF. This 
strategy requires much caution, as there can be differ¬ 
ences between the effects of the AO applied to the science 
data cube and to the standard star data cube. Neverthe¬ 
less, our tests revealed that, if such differences existed, 
in this case, they were of second order and, as a conse¬ 
quence, the Richardson-Lucy deconvolution not only did 
not compromise the data cube of M104, but actually im¬ 
proved significantly its spatial resolution. A comparison 
with an /-band image of M104 obtained with the Wide- 
field Planetary Camera 2 (WFPC2) of the HST revealed 
that the FWHA1 of the PSF of the data cube was ~ 0"3, 
before the Richardson-Lucy deconvolution, and ~ 0"15, 
after the Richardson-Lucy deconvolution. For more de¬ 
tail s about the da t a treatment procedure we performed, 
see IMenezes et alJ (12014120151) . An image of the treated 
data cube of M104 collapsed along the spectral axis, to¬ 
gether with its average spectrum, is shown in Figure [l] 

The main feature in the image of the collapsed data 
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Figure 1. Left: image of the treated data cube of M104 collapsed along the spectral axis. Right: average spectrum of the treated data 
cube of M104. 


cube of M104 is an elongated structure with a position 
angle (PA) of 63° ± 4°. This is the inner stellar disk 
of this galaxy, which was analyzed in prev ious studies 
(|Burkheadlll986tlEmsellem et al.l[T994l 119961) . 

3. ANALYSIS OF THE EMISSION LINE SPECTRUM 

One of our interests in the data cube of M104 is the 
analysis of the emission lines in the nuclear region of 
this galaxy. However, this analysis requires an accurate 
starlight subtraction. In order to perform this procedure, 
we applied the Pen a lized Pixel Fitting (pPXF) method 
(iCappellari & Emselleml l2004h to all the spectra of the 
data cube. This procedure uses a combination of tem¬ 
plate spectra from a given base convolved with a Gauss- 
Hermite expansion, to obtain the final synthetic spec¬ 
trum. In this case, we used a base of stellar spectra 
observed with NIFS, described in detail in IWinge et all 
( 2009 1. Since this base has the same spectral resolution 
of our data cube, it is very appropriate to be used to 
fit the absorption lines in the data cube. On the other 
hand, this base does not provide a good fit for the stellar 
continuum. Therefore, the algorithm we used to apply 
the pPXF added a fourth degree Legendre polynomial to 
the combination of the template spectra from the base, 
in order to fit the stellar continuum. 

Before applying the pPXF to the data cube of M104, 
we shifted the spectra to the rest-frame, using z = 
0.003416 (Nasa Extragalactic Database - NED), and 
sampled them with AA = 1 A per pixel. All these steps 
were performed using scripts written in IDL. The pa¬ 
rameters provided by the pPXF are: the stellar radial 
velocity (V - *), the stellar velocity dispersion (er*), and 
the Gauss-Hermite coefficients /13 and /14. The Gauss- 
Hermite coefficient /13 reveals asymmetries in the stel¬ 
lar absorption lines while the Gauss-Hermite coefficient 
/14 ind icates deviati ons from a Gaussian pro file in these 
lines dGerhardl 119931 Ivan der Marel fe Franxl I1993H . A 
synthetic stellar spectrum is also provided by the pPXF 
method. The starlight subtraction we required for a de¬ 
tailed analysis of the emission lines in the data cube of 
M104 was performed by subtracting the obtained syn¬ 
thetic stellar spectra from the observed ones. Figure [2] 
shows an example of the fit (together with the fit resid¬ 
uals) obtained with the pPXF for a spectrum from a 
circular area, with a radius of 0"084, centered on the 


nucleus of the stellar bulge. 

A simple analysis of Figure [ 2 ] reveals that the fit re¬ 
produced, with good precision, the continuum and the 
absorption lines. The fit residuals show only one signif¬ 
icant emission feature at 2.1218 /im, which corresponds 
to the H 2 A 212 I 8 emission line. We believe that other H 2 
lines, like H2A22234 or H2A22478, were not detected in 
this case because they are immersed in the noise of the 
spectrum. The H 2 A 212 I 8 emission line was detected here 
probably because it is usually the most intense H 2 line in 
this spectral region. We extracted the total spectrum of 
the data cube, after the starlight subtraction, and then 
we calculated the integrated flux of the H 2 A21218 line, 
obtaining a value of Fh 2 = ( 6.8 ± 0.7) x 10 ~ 15 erg cnr 2 
s -1 . Since we observed only the H 2 A 212 I 8 emission line 
in the data cube of M104, we could not perform a reli¬ 
able analysis involving the ratios of H 2 lines, in order to 
determine if this molecular emission is due to thermal or 
non-thermal processes in this galaxy (lMourilll994f) . 

In order to analyze the morphology of the H 2 A 212 I 8 
emitting regions, we constructed an image of this emis¬ 
sion line. We also constructed an RGB composite image 
(based on the radial velocity values) for this line. The 
results are shown in Figure [ 2 j We can see that most 
of the H 2 A 212 I 8 emission comes from a structure with 
PA = 65° ± 9°, which is compatible (at la level) with 
the PA of the inner stellar disk. This, together with the 
observed morphologies, suggest that the molecular gas 
is disposed along the inner stellar disk. The center of 
this molecular gas distribution is located at a projected 
distance of ~ 0"12 (toward east) from the stellar bulge 
center. Figure [2] also shows an amplitude/noise (A/N) 
map of the H 2 A 212 I 8 emission line. We can see that the 
A/N values are in the range of 10 - 14 near the stellar 
bulge center and decrease significantly at farther areas. 
They are smaller than 3 at the most peripheral regions. 
In order to construct a velocity map for the H 2 molec¬ 
ular gas, we fitted Gaussian functions to the H 2 A21218 
emission line. However, with these low A/N ratios, the 
Gaussian fits did not provide reliable values for the ra¬ 
dial velocity. We performed a spatial binning of the data 
cube but, in order to obtain sufficiently high A/N ratios 
for the H 2 A 212 I 8 emission line, the spatial resolution and 
the number of resolution elements were severely reduced. 
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Figure 2. (a) Spectrum from a circular area, with a radius of 0"084, centered on the nucleus of the stellar bulge in the data cube of M104. 
The fit provided by the pPXF is shown in red and the fit residuals are shown in green, (b) Magnification of the same fit residuals shown 
at left, (c) Image of the H2A21218 emission line of the data cube of M104, after the starlight subtraction, (d) RGB composite image of 
the H2A212I8 emission line. The colors blue, green, and red correspond to the velocity ranges —438 < V r < —170, —156 < V r < 156, and 
170 < V r < 438 km s 1 , respectively, (e) Amplitude/noise ratio map of the H2A212I8 emission line. 


Considering all of that, we decided to not include the ve¬ 
locity map of the H 2 A 212 I 8 emission line in this paper. 

Although we could not construct a reliable kinematic 
map of the H 2 A 212 I 8 emission line, the RGB compos¬ 
ite image in Figure [2] indicates that the kinematics of 
the molecular gas seems to be compatible with a rota¬ 
tion around the nucleus. Therefore, we conclude that 
the obtained results indicate the existence of a rotating 
torus/disk structure, approximately co-planar with the 
inner stellar disk around the nucleus of Ml04. 

4. ANALYSIS OF THE STELLAR KINEMATICS 


Our second interest in the data cube of M104 is the 
analysis of the stellar kinematics around the nucleus. As 
mentioned before, the pPXF method provides the values 
of Vt, (7 *, /13, and /14. Therefore, since this method was 
applied to each spectrum of the data cube of M104, we 
constructed maps of these parameters, which are shown 
in Figure[3] In addition, in order to evaluate if the values 
obtained for the kinematic parameters are reliable, we 
constructed an signal-to-noise ratio (S/N) map (for the 
wavelength range of 2.2730 - 2.2777 /irn). The result is 
also shown in Figure [3] It is important to mention that 
the V* map was calibrated in order to ensure that V* was 
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Figure 3. Kinematic maps obtained with the pPXF method applied to the data cube of M104. An S/N ratio map for the wavelength 
range of 2.2730 - 2.2777 fim is shown at the bottom. 


equal to 0 at the stellar bulge center. In other words, the 
values of F» shown in the map in Figure [3] are the stellar 
radial velocities relative to the stellar bulge center. 

An analysis of the V* map in Figure [3] reveals the exis¬ 
tence of a rotating stellar disk around the nucleus, with 
a maximum velocity of F* (max) = 256 ± 7 km s -1 and a 
minimum velocity of F*(min) = —223 ± 7 km s -1 . The 
cr* map is a little noisier than the F* map and shows a 
peak with a„(peak) = 272 ± 9 km s — , which is located 
at a projected distance of ~ 0"21 (toward east) from the 
center of the F* map. The cr* values are smaller along the 
stellar disk than in regions at larger distances. This is 
evidence that the superposition of a “cold” rotating disk 
(with lower cr* values) and a “hot” bulge (with higher cr* 
values) is a good description for the stellar kinematics 


in the nuclear region of M104. Although the h 3 map is 
noisier than the previous two, it is possible to detect a 
pattern, as the smallest values of /13 (with a minimum of 
h 3 (min) = —0.18 ± 0 . 02 ) are located in the region with 
the largest stellar velocities and the largest values of /13 
(with a maximum of hs(max ) = 0.14 ±0.04) are located 
in the region with the smallest stellar velocities. Finally, 
the /14 map does not show clear structures or patterns. 

Figure[3]also shows that the S/N values are larger than 
30 in most of the FOV of the data cube (they are even 
larger than 100 near the center of the stellar bulge). Only 
at the most peripheral areas the S/N values are smaller 
than 30. This indicates that the features shown in the 
kinematic maps (like the off-centered cr* (pea/c), for ex¬ 
ample) are not artefacts caused by small S/N values. In 
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Figure 4. Left: image of the symmetric component (with a Sersic index of 3.7) of the nuclear region of M104, obtained with the 
decomposition of an HST image, in the I band. Center: image of the asymmetric component of the nuclear region of M104, obtained 
with the decomposition of the same HST image. Right: original HST image that was decomposed in a symmetric and in an asymmetric 
component. In all images, the orientation is the same as in the data cube of M104. 


order to evaluate if some of these kinematic features were 
introduced by the Richardson-Lucy deconvolution or by 
the spatial re-sampling of the data cube of M104, we 
applied the pPXF method to the non-deconvolved data 
cube, with spaxels of 0"021, and also to the spatially 
binned non-deconvolved data cube, with spaxels of 0"1. 
The obtained maps had the same characteristics of the 
ones shown in Figure [3] Therefore, we conclude that 
the Richardson-Lucy deconvolution and the spatial re¬ 
sampling did not affect the measurement of the stellar 
kinematics. Finally, we also extracted the kinematics of 
the data cube of M104 using a Gaussian line of sight ve¬ 
locity distribution. Again, the main features of the V* 
and (j* maps (including the off-centered cr*(peafc)) were 
compatible with the ones in Figured] 

The asymmetry in the R* map, together with the fact 
the the a*(peak) is displaced from the center of the R* 
map, are evidences for the existence of a rotating excen- 
tric stellar disk in the nuclear region of M104. In order 
to obtain an estimate for the mass of the SBH in this 
galaxy, we performed a simple dynamical modeling, as¬ 
suming a thin eccentric stellar disk around the nucleus. 
For this procedure, first of all, we superposed 126 con¬ 
centric elliptic orbits, taking as free parameters the ar¬ 
gument of the pericenter uj, the mass of the central black 
hole M., the inclination of the disk i, and the eccentricity 
of the disk e. The longitude of the ascending node was 
measured in the R* map (fl = 56°4); therefore, we kept 
this parameter fixed in our model. The calculation of 
the velocities in each one of the orbits also required the 
knowledge of the stellar mass within each orbit. In order 
to determine that, we retrieved an image of the nuclear 
region of M104, obtained with WFPC2 of the HST, us¬ 
ing the F814W filter. This image was calibrated in order 
to provide the values of the luminosity in the I band. 
After that, we decomposed the image into a symmetric 
component (representing the central region of the stel¬ 
lar bulge of M104) and an asymmetric one (representing 
the rotating stellar disk). We verified that the best re¬ 
sult of this decomposition is obtained when we assume 
a Sersic profile, with a Sersic index of 3.7, to the central 
region of the stellar bulge. Figure 0] shows the HST im¬ 
age, in the F814W filter, of the nuclear region of M104, 
together with the images obtained, after the decomposi¬ 
tion, for the symmetric and asymmetric components. In 
order to de-project the symmetric component, we applied 


a Multi-Gaussian Expansion (Cappellari 2002) and, as a 
result, obtained the values of the luminosity density (in 
Liq/ pc 3 ). Considering that, we were able to determine 
the stellar mass within each elliptical orbit of our model 
by assuming a mass to light ratio in the I band (M/Lj), 
which was taken as another free parameter. Once all 
the concentric orbits were superposed, the velocity val¬ 
ues corresponding to each one of them were calculated 
and the orbits were projected on the plane of the sky. 
As a result, we obtained a velocity map. 

The next step in our dynamical modeling was the 
construction of two synthetic data cubes: one represent¬ 
ing the rotating stellar disk and the other representing 
the stellar bulge. In order to do that, we first obtained 
a representative spectrum for the data cube of M104, 
which was taken as a linear combination of the spectra 
of the base used before to apply the pPXF method. 
The coefficients of this linear combination were provided 
by the pPXF applied to the average spectrum of the 
data cube of M104. The representative spectrum was 
obtained with a spectral range of 2.25 - 2.40 /xm. The 
synthetic data cube representing the rotating stellar 
disk was constructed by shifting the representative 
spectrum using the values of radial velocities given by 
the synthetic velocity map obtained before and then 
by placing each shifted spectrum in the position of the 
corresponding spaxel of the synthetic velocity map. All 
these spectra were also convolved with an estimate of 
the stellar velocity dispersion of the rotating disk (ad), 
which is another free parameter of our model. It is 
important to emphasize that we assumed a constant 
stellar velocity dispersion for the rotating stellar disk. 
The synthetic data cube representing the stellar bulge 
was constructed by placing the representative spectrum, 
with a radial velocity equal to 0, in the position of 
each spaxel. All these spectra were convolved with an 
estimate of the stellar velocity dispersion of the bulge 
(ab), which is also a free parameter in our model. In 
other words, we assumed that the stellar bulge shows 
a radial velocity equal to 0 and also has a constant 
velocity dispersion. Using the images of the symmetric 
and asymmetric components of the nuclear region of 
M104, obtained from the decomposition of the HST 
image described before, we were able to determine the 
relative weights between the spaxels of the constructed 
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Figure 5. Left: curves of V*, cr*, and /13, extracted along the kinematic axis. The observed curves are the points with error bars and 
the simulated curves of the best obtained model are shown in red. Center: simulated maps of V *, cr*, and h% of the best obtained model. 
Right: observed maps of V*, cr*, and /13, with a reduced FOV. 


synthetic data cubes. With these relative weights, 
the synthetic data cubes representing the rotating 
stellar disk and the stellar bulge were combined into 
one. Finally, each image of the combined data cube 
was convolved with an estimate of the PSF of the 
observation. In order to obtain the simulated kinematic 
maps, we applied the pPXF to the synthetic data 
cube, completing the procedure. In summary, the free 
parameters of our model are: w, M., i, e, M/Lj, ad, 
and a;,. Since the /14 map shown in Figure [3] is quite 
noisy and does not show clear structures, we only tried 
to reproduce the maps of the values of 14 , cr*, and ft.3. 
Considering that we are more interested in the stellar 
kinematics near the central black hole and also in order 
to avoid higher imprecisions in the calculations of the 
stellar mass, only a region within a radius of 0"85 was 
simulated. The free parameters of our model were varied 
and this process was repeated several times, in order to 
minimize the value of the total y I 2 (x§-), which is given by 


Xt = 






I \_l_f J V 

'■Min( x £ 3 r 


(1) 


where Xyj xt ,, and Xh 3 are the x 2 s calculated for 
the simulated maps of 14 , cr*, and / 13 , respectively, 
and Min(xvJ, an( t Min(Xh 3 ) are the min¬ 

ima values obtained for \v, 1 Xa, 1 an( t Xh 3 > respectively. 
Our tests revealed that, in our dynamical modeling, the 
weighted mean given by equation ( 1 ) is the most effective 
way to combine the values of Xv > Xct, > an d Xh i n t° one - 
We defined the x 2 calculated for the simulated map of 
a kinematic parameter A * in the following way: 


2 _ 1 \~~\Ny Wij -Iij -(Aij (observed)—Aij (simulated)) 

Xa* — 7 2_,i=i 4*7=1 7^ 


(2) 
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Table 1 

Parameters of the Best Model Obtained to Reproduce the Stellar 
Kinematics of the Data Cube of M104. 


Parameter 

Value 

M. 

(9.0 ± 2.0) x 10 8 M 0 

e 

0.12 ±0.05 

i 

80° ± 2° 

CO 

5° ±9° 


(260 ± 14) km s _1 

°d 

(120 ±8) km s _1 

M/L / 

3.0 ±0.5 


where A^ (observed) and Aij (simulated) are the values 
of the spaxel (i, j) of the observed and simulated maps 
of A *, respectively, N x and N y are the number of spaxels 
along the horizontal and vertical axis, respectively, 
is the uncertainty of A * in the spaxel Wij is a step 
function equal to 1 for areas closer than 0"20 to the kine- 
matical axis and equal to 0 for areas at larger distances, 
Iij is the integrated flux of the data cube of M104 in 
the spaxel (i,j), and / is the total integrated flux of the 
data cube of M104. The values of the uncertainties aA tj 
of the kinematic parameters analysed in this work (V*, 
cr*, and / 13 ) were obtained via Monte Carlo simulations. 
The calculations of the values of Xa involved the inte¬ 
grated fluxes of the data cube of M104 (J^ and I) and 
also the step function wij, in order to give more weight 
to the areas closer to the kinematical axis of the data 
cube. If we had not adopted this strategy, most of the 
weight in the calculation of x'a would have come from 
the outer parts of the image (with a higher number of 
spaxels). However, we are more interested in reproduc¬ 
ing the kinematic parameters in areas near the kinematic 
axis. Besides that, since we are applying a model of a 
thin disk, it is expected that the limitations of our model 
be more evident in areas farther away from the plane of 
the disk. Therefore, considering that the rotating disk 
of M104 is almost edge-on, the limitations of our model 
would be more significant in regions at larger distances 
from the kinematic axis, originating higher discrepancies. 
Based on all of that, we decided to perform the calcula¬ 
tion of the values of Xa giving more weight to areas 
near the kinematic axis. There are different ways to give 
more weight to the central regions of the images during 
the calculation of the y 2 ; however, once again, our tests 
revealed that the specific calculation given by equation 
( 2 ) provides the best results in this case. 

The simulated maps, and the simulated curves ex¬ 
tracted along the kinematic axis, of 14 , cr*, and /13 of 
the best obtained model, with Xt = 1-519, are shown in 
Figure 0 In the same Figure, we show the correspond¬ 
ing observed maps (with reduced FOV). The kinematic 
parameters of the best model are shown in Table Q] We 
estimated the uncertainties (ltr) of th e para meters by fol- 
lowing the approach described in Ide Rosa et all (120Ilf) : 
we plotted an histogram for each parameter of the simu¬ 
lation, considering only cases with Xt — Xmm < 1- Then, 
a Gaussian function was fitted to each histogram. The 
value obtained for the square deviation of each Gaus¬ 
sian was taken as the uncertainty of the corresponding 
parameter. 

An analysis of Figure 0 reveals a good agreement be¬ 
tween the observed and simulated 14 curves, as all the 


simulated values are compatible with the observed ones 
at lcr or 2a levels. There is also a good agreement be¬ 
tween the simulated and observed 14 maps within a dis¬ 
tance of ~ 0"2 from the kinematic axis. On the other 
hand, we can see a significant discrepance between the 
simulated and observed values in areas at distances of 
0"2 — 0"4 from the kinematic axis. As mentioned before, 
this is probably a consequence of the limitations of a thin 
disk model. For regions at larger distances, both the sim¬ 
ulated and observed maps show values close to 0. Almost 
all the observed values of the cr* curve are compatible, at 
lcr level, with the simulated ones. The general morphol¬ 
ogy of the cr* map, including the position of the peak, is 
well described the model. However, again, we can see cer¬ 
tain discrepancies, most likely related to the limitations 
of the model, in areas at distances of 0"2 — 0"4 from the 
kinematic axis. For regions at larger distances, the ob¬ 
served cr* map is much more irregular. It is probable that 
these irregularities represent deviations from the model 
we used here and that could explain the observed discrep¬ 
ancies. In the case of the values of ft. 3 , the situation is 
somewhat different. Although the general patterns of the 
observed /13 curve and /13 map are reasonably reproduced 
by the model, there are many discrepancies that should 
be discussed. The points with positive values of the sim¬ 
ulated and observed /13 curves are compatible at lcr or 2 cr 
levels. However, most of the other points of these curves 
are not compatible, even at 3 a level. The observed /13 
map is quite noisy and shows several assimetries and ir¬ 
regularities that were not adequately reproduced by our 
model. Part of these irregularities may be the result of 
uncertainties of the pPXF method. However, other irreg¬ 
ularities may be real and, therefore, represent deviations 
from our thin disk model. One example is the central 
features in the h .3 map correlated to the 14 map. As 
explained above, the S/N values in the areas near the 
stellar bulge center are significantly large. In addition, 
the /13 maps obtained with the pPXF applied to the non- 
deconvolved data cube, with spaxels of 0 " 021 , and to the 
spatially binned non-deconvolved data cube, with spax¬ 
els of 0 " 1 , also shows these structures correlated to the 
14 map. This indicates that these features are not arte¬ 
facts introduced by the Richardson-Lucy deconvolution 
or by the spatial re-sampling. Based on this discussion, 
we conclude that the central structures in the /13 map 
are real and not a consequence of uncertainties caused 
by small S/N values or an artefact introduced by the 
methodologies used during the data treatment. These 
/13 features correlated to the 14 map are not expected 
for a rotating disk within a bulge, which suggests the 
existence of deviations from the model we applied here. 
The high asymmetry in the h3 map is probably also a de¬ 
viation from the model, as this behavior was observed in 
all /13 maps obtained with the applications of the pPXF 
to the non-deconvolved data cubes (with and without a 
spatial binning) of Ml04. 

5. DISCUSSION AND COMPARISON WITH PREVIOUS 
STUDIES 

In iMenezes et al.l (|2013| ). we have detected a dusty 
dense structure (with PA = 72° ± 14° and an outer radius 
of 0"40, corresponding to r = 18 pc) in the nuclear region 
of M104, which we interpreted as an obscuring torus/disk 
that collimates the AGN emission of this galaxy. The PA 
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of this torus/disk is compatible, at ler level, with the PA 
of the rotating H 2 torus/disk (PA = 65° ± 9°) observed 
in this woA. This sugge sts that the structures detected 
in iMenezes et al.1 ( 2013 1 and in this work are coplanar 
and associated with the same obscuring torus/disk. 

The morphology of the 14 map we obtained for the nu¬ 
clear region of Ml04 is compatible with the m orphology 
of the map obtained bv lEmsellem et all (Il996tl . using the 
integral-field spectrograph TIGER. The morphologies of 
our < 7 * and /13 maps are also similar to what was obtained 
by these authors, but some discrepancies can be seen. We 
believe that these discrepancies can be explained by the 
fact that the spatial resolution of our NIFS data cube 
(~ 0"15) is significantly higher than th e spat ial reso lu¬ 
tion o f the TIGER data cube obtained by Emscllc m et al.1 
(il996l) (~ 0"9). Therefore, we are able to detect struc¬ 
tures and irregularities that are impossible to be observed 
in TIGER data. 

The observed rotation curve shown in Figure [5] is 
consis tent wit h the curves obta i ned b y [Emscllc m et al.l 
(11996f l and bv iKormendv et all (1199611 . These previous 
works also obtained cr* curves, which are similar to our 
cr* curve, but some differences can be easily detected. 
Our cr*(peafc) (= 272 ± 9 km s -1 ) is nearly compati- 
ble, at 3 a level, with the a^(peak) of lEmsellem et all 
(fT996tl (= 248 km s 1 ) and is not compatib le with the 
cr*(peafc) obtained bv IKormendv et all (|1996l ) (= 321 ±7 
km s _1 ). We believe that these discrepancies are caused 
by differences in the spatial resolutions of these works. 
One interest ing topic is that t he po sitio ns of a*(peak) de¬ 
termin ed bv lEmsellem et akl (1199611 and IKormendv et al.1 
(1996) are apparently coincident with the position of the 
14 map center. However, in the present work, we ob¬ 
served a projected distance of ~ 0"21 between the posi¬ 
tion of er*(peafc) and the 14 map center. T he separatio n 
of the se points could not be detected bv lEmsellem et all 
(Il996ll due to the spatial resolution of the observation. 
On the other hand, the non-detection of this behavior 
of < 7 * (peak) by IKormendv et al.l (11996H is more difficult 
to explai n. Three points of the cr* curve of Ml04 pre¬ 
sented bv IKormendv et akl ( 1996 1 were obtained with the 
Faint Object Spectrograph (FOS) of the HST (while the 
other points were obtained with the Subarcsecond Imag¬ 
ing Spectrograph of the CFHT). These three FOS values 
of cr* were determined from spectra resulting from three 
observations, using the 0"21 square aperture, at distances 
of 0"032 W, 0"163 W, and 0"207 E from the nucleus. All 
these observations were taken along the east-west direc¬ 
tion. Based on this information, we conclude that the 
spectrum of the region at 0"207 E from the nucleus could 
have revealed the off-centered cr*(peafc) we observed in 
the data cube of Ml04. Howeve r one should note that 
the cr* curve in IKormendv et al.l ( 1996 ) reveals that the 
difference between the < 7 * values provided by the spectra 
taken at 0'.'032 W and at 0"207 E from the nucleus are 
~ 10 km s -1 . On the other hand, the error bars of the 
point at a radius of 0"032 W are ~ 25 km s _1 while the 
error bars of the point at a radius of 0"207 E are ~ 10 
km s -1 . Therefore, considering the difference between 
the cr* values in these two points and the error bars, we 
conclude that the separation between the <7*(peafc) and 
the stellar bulge center we observed in the data cube of 
Ml04 is actually not incompatible with the results ob¬ 


tained by IKormendv et al.l (|1996T) . There are also sig¬ 
nificant di screp anc ies betw een the /13 curves obtained by 
us and bv lEmsellem et akl (|1996f ). However, considering, 
again, the differences in the spatial resolutions of the cor¬ 
responding observations, we can say that the observed 
discrepancies are not surprising. The central features in 
the /13 map correlated to the 14 map may be i n dicative 
of a bar in the central region of Ml04 dChung fe Bureaul 
12004 llannuzzi fe Athanassoulal I2015D . The hypothesis 
of t he existence o f a ba r in Ml04 was already discussed 
by lEmsellem & Ferruiil ( 2000 ! and should be taken into 
account in future studies. 

The mass of the SBH of M104 we obtained with a sim¬ 
ple model of a thin eccentric disk (9.0 ± 2.0 x 10 8 Mq) 
is compat ible, at ler l evel, with the value of 10 9 Mg esti¬ 
mated bvlKormendvl (11988ft . lEmsellem et al.1 (|1994fi . and 
IKormendv et al.l (119961) . Our result is also compatible, 
at 2cr level, with the masses o btained b y Jarde Fet al.1 
(120111) (6.6±0.4x 1O 8 M 0 ) and bv lMagorrian et al.1 (I1998T) 
(6.47±g;?| x 10 s M 0 ). The value of M/Lj resulting from 
our mode l is compatible, at ler level, with the value ob¬ 
tained by IJardiT^al] (j201l]) (M/Lj = 3.40±0.05). Our 
model is certainly much more simplistic than other dy¬ 
namical models, like the Scliwarzschild and the Jeans 
methods, which reproduce in detail the gravitational po¬ 
tential around the SBH. The Schwarzschild method in¬ 
volves the determination of the gravitational potential 
using the brightness profile of the galaxy, the construc¬ 
tion of an orbit library, the projection of the orbits on the 
plane of the sky, and then the superposition of the orbits, 
in order to reproduce the observed kinematics. Usually, 
the Schwarzschild method is applied assuming that the 
galaxy is axisymmetric or triaxial. The Jeans method 
usually assumes axisymmetry and involves the solution of 
the Jeans equations, in order to obtain the gravitational 
potential and the kinematics around the SBH. Although 
our method is not as detailed as these two, it has certain 
advantages. The Schwarzschild method, for example, is 
frequently affected by degeneracies and difficulties in the 
convergence. Such problems are rarely observed in the 
model of a thin eccentric disk, due to its simplicity. The 
Jeans method often assumes axisymmetry, which is not 
required by the thin disk method. We have already used 
t he implementat ion of the Jeans method developed by 
Cappellarij (120081) to try to simulate the k inematic maps 
of the data cube of M104 (|Menezesll2012D . However, we 
were not able to reproduce certain aspects of these maps, 
like the displacement of <r*(peak) from the U* map center. 
This can be explained by the fact that the implementa¬ 
tion of the Jeans method we used assumes axisymmetry, 
which is not the most appropriate approximation to be 
used for this galaxy. As we verified in this work, the as¬ 
sumption of a rotating eccentric disk is more adequate. 
In summary, the disadvantages of our model of a thin 
excentric disk are: it, obviously, assumes the existence 
of a thin rotating disk, which is not necessarily the case; 
it is less detailed than the Schwarzschild and the Jeans 
methods. On the other hand, its advantages are: due 
to its simplicity, it is usually not affected by degenera¬ 
cies and problems with convergence; it does not require 
axisymmetry. The fact that the mass of the SBH we ob¬ 
tained is compatible, at ler or 2 er levels, with the values 
found in previous studies using very detailed dynamical 
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models indicates that the model of a thin eccentric disk 
can be used to describe, with good precision, the stellar 
kinematics around the SBH in M104. 

The presence of an eccentric disk around a SBH of a 
galaxy is not actually new, as this structure was already 
detected in the nuclear region of M31, the Andromeda 
galaxy. The eccentric disk in this object is responsi¬ 
ble for the existence of a double nucleus (lLauer et al.1 
119931) . The two components are called PI (the bright¬ 
est one) and P2 (located, approximatelly, at the cen¬ 
ter of the stellar bulge). PI coincides with the apocen- 
ter of the stellar dis k, while P 2 is located in the vicin¬ 
ity of the SBH. [Bender et al.l (|2005H verified that the 
SBH is actually located in a structure embedded in P2, 
called_P3j which is probably a cluster of A type stars. 
lLauer et al.l ( 2012 ) confirmed that P3 is a cluster of blue 
stars around the central black hole. Different studies de¬ 
tected an asymmetry in the stellar rotation curve of the 
nuclear region of M31 and also a displacement between 
the kinema tic center and the peak of the stellar velocit y 
dispersion dBacon et al.lll994120011 : [Bender et al.ll2005ll . 
This last feature is located in an area close to the peri- 
center of the disk, which shows a lower surface bright¬ 
ness than the area corresponding to the apocenter (PI). 
All these photometric and kinematic characteristics were 
successfully reproduced with models involvin g a rotating 
eccentric stellar disk around the SBH rtTremainel 119951 : 
iSalow fc Stalled 12004 iPeiris fc Tremaine! 120031) . In the 
case of M104, our dynamical modeling revealed that the 
eccentricity of the stellar disk is smaller than the val¬ 
ues obtained for M31 (which can be higher than 0.3 for 
certain parts of the disk, as revealed by the studies men¬ 
tioned before). Nevertheless, we can identify many sim¬ 
ilarities in the stellar kinematics and in the photometry 
of these two objects. The kinematic similarities are the 
asymmetry in the rotation curve and the displacement 
between the kinematic center and the peak of the stellar 
velocity dispersion. In the case of the photometry, the 
similarities are more difficult to detect. A careful analysis 
of the HST images in Figure [4] reveals a slight decrease 
in the surface brightness of the disk in an area at a dis¬ 
tance of ~ 0"25 eastward from the nucleus. This region 
corresponds approximatelly to the position of cr*(peak) 
and, according to our dynamical modeling, it also co¬ 
incides with the pericenter of the disk. On the other 
hand, the surface brightness of the disk is higher in the 
diametrically opposed area, which, according to our dy¬ 
namical modeling, coincides with the apocenter of the 
disk. Since the purpose of this work is not to analyze 
in detail the photometry of the nuclear region of M104, 
we will not give further detail about this topic. How¬ 
ever, the discussion above reveals that the photometry 
and the stellar kinematics in the nuclear regions of Ml04 
and M31 (which shows a well studied rotating eccentric 
stellar disk) have certain similarities, although some of 
them are subtle. This indicates that the model of an 
eccentric stellar disk is actually a very natural choice to 
reproduce the stellar kinematics in the nuclear region of 
M104. 

6. SUMMARY AND CONCLUSIONS 

We analyzed a data cube of the nuclear region of M104, 
the Sombrero galaxy, obtained with NIFS on the Gemini- 
north telescope. After a subtraction of the stellar con¬ 


tinuum, we detected the H 2 A 21218 emission line in the 
spectra of the data cube. The image, and also the corre¬ 
sponding RGB image, of this line indicate that the molec¬ 
ular gas is disposed along a rotating torus/disk structure, 
aproximately co-planar with the inner stellar disk around 
the nucleus of Ml04. This molecular structure is appar¬ 
ently coplanar with anoth er s t ructu re, made of dust, de¬ 
tected by us in lMenezes et al.l (12013D . which we identified 
as an obscuring torus/disk. Therefore, we conclude that 
these two findings are probably associated with the same 
obscuring torus/disk. 

Using the pPXF method in the data cube of M104, we 
obtained maps of U*, <r*, / 13 , and / 14 . We detected a ro¬ 
tating stellar disk around the nucleus. We also observed 
that the < 7 * values along the stellar disk are consider¬ 
ably smaller than the values at the surrounding areas. 
This indicates that the stellar kinematics in the nuclear 
region of M104 can be described by the superposition 
of a “cold” rotating disk (with smaller <r* values) and a 
“hot” bulge (with larger cr* values). Assuming the ex¬ 
istence of a thin eccentric disk, we applied a dynamical 
modeling, in order to reproduce the V*, cr*, and h 3 maps 
and also obtain an estimate of several kinematic param¬ 
eters. The value of the mass of the SBH we obtained 
(M. = 9.0 ±2.0 x 10 8 Mq) is compatible, at lcr or 2 a lev- 
els, with the estimates found by many previous studies 
(lKormendvlll988tlEmsellem et al.lll994lKormendv et al.1 

119901 : iMagorrian et, al.l 119981: I.Iardel et al.l 1201 111 . Our 
modeling reproduced the main properties of the V* and 
cr* maps, specially within a distance of 0"2 from the kine¬ 
matic axis, and also the general behavior of the /13 map 
(although the discrepancies, in this case, are consider¬ 
ably larger). For distances of 0"2 — 0"4 from the kine¬ 
matic axis, the discrepancies between the simulated and 
observed kinematic maps increase substantially, which is 
probably a consequence of the limitation of a thin eccen¬ 
tric disk model. Nevertheless, considering all the results 
and also the compatibility between the value of M, ob¬ 
tained by us and by many previous studies, we conclude 
that the model of thin eccentric disk is a reasonably good 
description for the stellar kinematics around the SBH in 
M104. 

The model of an eccentric disk around the SBH of a 
galaxy was already applied to reproduce the stellar kine¬ 
matics around the nucleus of M31. Considering that, 
and also the results obtained in this work, we believe 
that this model may also be used to reproduce the stel¬ 
lar kinematics in the nuclear region of other galaxies and 
obtain reliable estimates for the masses of their central 
SBHs. 


Based on observations obtained at the Gemini Ob¬ 
servatory (processed using the Gemini IRAF package), 
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tional Research Council (Canada), CONICYT (Chile), 
the Australian Research Council (Australia), Ministerio 
da Ciencia, Tecnologia e Inovagao (Brazil) and Ministe¬ 
rio de Ciencia, Tecnologia e Innovation Productiva (Ar¬ 
gentina). We thank FAPESP for support under grants 
2012/02268-8 (RBM) and 2011/51680-6 (JES) and also 
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